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Abstract 

Cold atoms are traditionally trapped by the dipole force in periodically spaced potential wells 
induced by the standing laser field. We derive here a theory beyond the conventional dipole 
approximation which provides field / atom coupling potential terms that so far have not been taken 
into consideration in theoretical or experimental studies. We show that for some atoms for specific 
laser parameters despite the absence of dipole transition laser trapping is still possible due to the 
quadrupole force. Illustrative numerical calculations for Ca and Na trapped by the quadrupole 
force are given. 



1 




lattices, has opened a new subject of research. For reviews 
^| and references therein. Laser cooling and the formation 



The possibility of trapping cold neutral atoms in periodic potentials formed by standing 
laser waves, so called optica 
see for example Refs. 

of optical lattices were explained by the ability to influence the translational motion of 
neutral atoms by the laser induced dipole force J]. Until now, the formation of optical 
lattices has been studied theoretically and experimentally exclusively due to the dipole 
interactions. An important question arises, however, on what happens when the dipole 
transition matrix element vanishes due to a symmetry property of the studied atomic system. 
This is the case, for example, when the laser frequency is detuned from an excited state which 
is associated with a d-type symmetry orbital whereas in the ground state the valence shell 
electron occupies an s-type symmetry orbital. Does this imply that in such a case atoms 
can not be trapped because one can not form an optical lattice ? 

The purpose of the present study is to demonstrate that it is possible to trap atoms 
even when the dipole transition between the two relevant electronic levels coupled by the 
laser is prohibited by symmetry. We derive here the Hamiltonian of optical lattices without 
imposing the dipole approximation, and study the possibility of generating optical lattices 
due to the quadrupole contributions to the trapping potential which until now have been 
neglected. Two illustrative examples are presented explicitly. Namely, we show that by 
using realistic laser parameters the quadrupole force induced optical lattice can be formed 
and trap Ca and Na atoms. 

Let us consider an atom interacting with a linearly polarized laser light. The light propa- 
gation direction is chosen to be x, and the field is assumed to oscillate along the z-direction. 
For the sake of simplicity and without loss of generality we describe the field free atomic 
Hamiltonian by an effective one electron model Hamiltonian, . The formulation given 
here can be extended in a straightforward manner to the general case of a many electron 
atom (or ion) interacting with an arbitrary (not necessarily linearly polarized) light pulse. 
Such a generalization will be elaborated elsewhere By switching from the laboratory 
frame coordinates to the relative q = (x,y,z) and center-of-mass (cm.) R = (X,Y,Z) 
coordinates, one obtains the full Hamiltonian of an effective one electron atom in a linearly 
polarized laser field, 

H(t) = ~ V| + H c 7(g) + V coup (q, R, t) (1) 
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with an interaction term \ coup (q, R,t) = e(m e c) _1 A m j> z + e 2 (2m e c 2 )~ 1 A$ + 
e(cM)" 1 AmPz+^ 2 {^{M —m e )c 2 )~ l Aj$ . Here, the A-factors are defined in terms of the time, 
relative coordinate, and cm. coordinate dependent vector potential A(x, X, t), as follows: 
A me = [A(xCi+X,t)+(C 2 /Ci)A(-xC2+X,t)} ; A M = [A(xZi+X,t)-A{-x£ 2 +X,t)] ; A® = 
A 2 {x^\ + X, t) and A M ^ = A 2 (—x^2 + X, t). The symbols m e and M stand, respectively, for 
the mass of the electron and the atom, and parameters £1 = (M — m e )/M , £2 = m e /M . 
When the field intensity I is weak (/ is proportional to the squared field amplitude) the lin- 
ear interaction terms A me (x, X, t)p z and Am(x, X, t)Pz become dominant over the quadratic 
interaction terms A$ and A$ so that the latter two contributions can be neglected. More- 
over, since M 3> m e it is clear that £1 ~ 1 and thus the dominant term in \ coup (q, R,t) 
is the first one. In summary, we conclude that the coupling between the cm. and relative 
coordinates can be considered as 

Vcou P (q,R,t) = e(m e c)~ 1 A me (x,X,t)p z . (2) 

In standard treatments of the atom-field Hamiltonian the dependence of the vector po- 
tential on the relative coordinate has been neglected [jj. That is, A(x +X,t) ~ A(X,t). 
Here, we do not use this so called dipole approximation. Instead, the spatial dependence 
of the vector potential is taken into account, via the Taylor expansion around x = 0. An 
important question arises on what should be taken as the small parameter in such a Taylor 
series expansion. In the case when the laser intensity is low, the associated time depen- 
dent electronic wavefunctions are only slightly different from their field-free counterparts. 
Therefore, they possess exponentially small values for q > ao, where ao stands for the size 
of the atom. For weak fields, as in our case, ao is essentially equal to the radius of the field 
free atom, tq. Of course, in the presence of very strong fields ao varies with the intensity 
as the electron oscillates with the field. Roughly speaking we may think that the relevant 
radius of the atom is then r + a where a = (eA ) 2 / (c 2 m e u L ) and A is the field amplitude 
^]. In order to find a qualitative criterion for convergence of our Taylor series expansion 
of A(x + X, t), the above estimated atomic radius ao should be compared with the wave- 
length A of the laser. Clearly, A should get sufficiently large values such that A ^> ao • The 
small parameter of the considered Taylor series expansion is thus found to be fciao, where 
k L = 2tt/A. 

Following the above discussion, when the wavelengths which constitute the dominant part 
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of the laser pulse are sufficiently large, A > ao, the two leading order terms in the Taylor 
series expansion of the vector potential are given by 

A(x + X,t) ~ A(X,t) + xdA{X,t)/dX. (3) 

The first term in Eq. Q is the dipole term, whereas the second one is the quadrupole term 
which is taken here into consideration for the first time in contrast to other theoretical 
studies published up to now. Substituting Eq. Q into formula (J2J) one finds that 

V coup (q, R, t) = V dip (z, X, t) + V qd (z, x, X, t) (4) 

where the dipole and quadrupole potential terms that couple the internal (electronic) and 
cm. degrees of freedom are respectively given by, ~Vdi p (z, X, t) = (e/c)m~ 1 A(X,t)p z and 
V qd {z,x : X,t) = {e/c)m- l [dA(X,t)/dX]xP z . 

Since the field oscillates periodically in time the Hamiltonian is time periodic with the 
period T = 2tt/ul where ujl is the laser frequency. In such a case the solutions of the time de- 
pendent Schrodinger equation are the Floquet solutions, similarly as the Bloch states are the 
solutions of the Schrodinger equation for spatially periodic potentials. Within the framework 
of the Floquet theory, the time averaged optical potential, V op t(X), is associated with the 
quasi-energy eigenvalue of the Floquet operator, Tip = —ihd/dt + H^ F (g) + \ coup (q, R,t), 
(for Floquet formalism see for example Ref.j^l), 

H F ^(q, X, t) = VoptiX^iq, X, t) . (5) 

The eigenvalue solution of Eq. (JHJ provides the definition of the optical lattice potential. 
Strictly speaking, optical potentials for the atoms exist only in the adiabatic approximation 
where one can separate the atomic center of mass motion from the electronic motion (simi- 
larly to the Born-Oppenheimer or adiabatic approximation which separates the nuclear and 
electronic motions in molecules). Within the adiabatic approximation, our claim that the 
optical potential equals to an eigen-energy of the Floquet operator is exact. Of course, there 
are infinitely many eigen-energies and hence optical potentials (similarly to molecules where 
in each electronic state there is a potential for the vibrations). In equation (jHJ) we have 
chosen that eigen-energy which correlates to the field-free atomic ground state. In the weak 
field limit we will show below that in the dipole case the same expression as that used in 
the literature results, whereas in the quadrupole case a novel expression for the quadrupole 
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force induced optical lattice potential is obtained. Note however that for sufficiently strong 
fields, the literature expression (which is found to be a second order perturbation theory 
result) fails whereas the eigen-energy of the Floquet operator still provides an exact optical 
potential. 

Let us derive now an analytical expression for the optical lattice potential which in- 
cludes the leading order term beyond the dipole approximation. From Eq. (J3J) it is 
clear that the only information we need in order to go beyond the dipole approxima- 
tion is A(X, t) and its first order derivative with respect to the cm. coordinate X. For 
atoms in a standing laser beam A(X,t) = A$ cos(fciX) cos(a;x,i), where ujl is the cor- 
responding laser frequency, and k L = u L /c. Suppose that the atomic ground state 
|g) (with energy £®) is significantly coupled by the laser light only to a single excited 
electronic state |e) (with energy £® = £® + 7kJ a tom)- By following the "traditional" 
approach based upon the adiabatic elimination scheme jjj], we obtain the time depen- 
dent optical lattice potential (including the new term beyond the dipole approximation), 
V opt (X,t) = \(e\y dip (z,X,t) + V qd (z,x,X,t)\g)\ 2 /(hA L ), where, A L = u atom - The 
characteristic timescale of the center of mass dynamics is by several orders of magnitude 
smaller than the time period of an optical cycle, T = A/c = lixjio^. Therefore, the time 
averaged optical lattice potential reads as V opt (X) = (1/T) Jq V opt (X, t) dt. In case that the 
laser detuning with respect to the atomic frequency is such that it couples an s-type atomic 
state |g) with a p-type atomic state |e), the quadrupole contribution to V op t(X) vanishes, 
and the optical lattice potential is reduced to the expression 

Vg(X) = /i|fi dip | 2 (4A L )^ cos 2 (A; L X) (6) 

where the Rabi frequency is given by Q dip = (e A ) / (hcm e ) (e|p z |g). This is, of course, the 
well known conventional expression for the optical lattice potential formed by the dipole 
force jij. 

Let us consider now the case when the laser detuning is off resonance but not from a 
p-type atomic state but from a c?-type atomic state. Consequently, the force exerted on the 
atom is of quadrupole nature since the dipole contribution vanishes. Note that for Calcium, 
for example, the lowest excitation is from an s to a d orbital. Introducing the "generalized" 
Rabi frequency Q q d = (e A )/(hcm e ) k L (e\xp z \g) we obtain the formula 

(4 Ax) -1 sin 2 (k L X) (7) 
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V f t (X) = h 



n 



qd 



which exhibits remarkable similarity to the standard expression (JBJ) for the dipole optical 
lattice potential. Before proceeding further, let us recall once again that both relations (JHJ) 
and (J7J) can be derived directly by applying the second order perturbation theory to the 
Floquet eigenproblem ©j^l- 

Our theoretical analysis concludes with a very brief discussion of the as yet neglected 
incoherent effects (friction force, spontaneous emission) which are caused by the quantized 
nature of the electromagnetic field. It is known that the mentioned dissipation processes 
can be safely ignored when ^> T, where T is the natural linewidth of the excited electronic 
state. For the dipole |e) — > |g) transitions, 

r dip = 4e 2 w atom (3m^ c 3 )- 1 |(e|p z |g)| 2 , (8) 

see, e.g., Ref. The derivation leading to Eq. (|SJ) can be extended in a straightforward 
manner to cover also the situation when the transition |e) — > |g) occurs solely via the 
quadrupole coupling. The final result is then given by formula 

T qd = 4e 2 ^ tom (5m^ C 5 )- 1 |(e|x^|g)| 2 . (9) 

Let us apply now our above theoretical results and study the feasibility of trapping 
Na and Ca atoms by the quadrupole force in the periodically spaced potential wells 
induced by the standing laser waves. Our model Hamiltonian consists of the ground S 
state, \Eg), and the lowest excited P and D states, \Ep) and \Ep>). The correspond- 
ing energy levels for Ca and Na were taken from the NIST atomic database. The 
field free electronic Hamiltonian is a diagonal 3x3 matrix, E e , where the diagonal 
terms are the mentioned atomic energies E$, Ep, Ep,. Following Eq. (jSJ), the interac- 
tion potential induced by the linearly polarized standing laser field is described by the 
symmetric 3x3 matrix \ coup (X,t), where the diagonal matrix elements vanish. The 
off diagonal matrix elements are connected either with the dipole or the quadrupole 
transitions, and take an explicit form: {V coup )i j 2 = eA /(m e c)(Es\p z \Ep) cos(/clX) cos^upt), 
(y coup ) 2 ,3 = eA /(m e c)(E P \p z \E D ) cos(k L X) cos(u L t) and {V coup )i ;i = 
— eA / (m e c)kL(E s \xp z \E D ) sm^k^X) cos(a>£t). The corresponding dipole and quadrupole 
transition matrix elements for Na and Ca have been calculated using equations (jHI) and 
(jHJ) from available theoretical/experimental data for the associated transition rates. More 
specifically, the dipole transition linewidths and T^p are extracted from the NIST 

6 



database, while the quadrupole terms T^f and T^f are taken from the literature Q|. 
Referring to previous discussion of Eq. (jSJ), we recall that there is a single quasi-energy 
eigenvalue of the Floquet hamiltonian Tip = E e + V coup (X, t) —ih(d/dt) which is associated 
with the field free ground state and which is interpreted physically as the appropriate 
numerically exact optical lattice potential. 

In the first step of our calculations we computed the variation of the depth of the optical 
lattice potential well V ma , x = Max[V op t(X)] for Na atoms as a function of the laser intensity, 
/ [W/cm 2 ] oc (kT,Aftl2) 2 for the laser parameters as were used recently in the experiment 
of Phillips group |l4j. That is, Na atoms are trapped in a one dimensional optical lattice 
induced by a laser with A = 589 nm, and the detuning from the first excited P state of Na is 
equal to A^ = 14 GHz. Our calculations show that for / < 100 W/cm 2 the analytical results 
derived from the second order perturbation theory, Eq. ©, are in a remarkable agreement 
with the numerically exact values. Therefore, indeed within this region of field intensities 
the optical lattice is formed exclusively by the dipole force. 

The situation is entirely different when the laser wavelength is taken to be A = 343 nm 
and the detuning from the first excited D state is Ai = 1 GHz. From Fig. 1 one can 
see that the quadrupole force induced optical lattice potential predicted by Eq. ((7|) is now 
dominant and in complete agreement with the numerically exact solutions of Eq|3 for the 
entire range of displayed field intensities. The dipole force contribution to the formation of 
Na optical lattice as calculated from Eq. © is, for the considered laser parameters, by one 
order of magnitude smaller than the contribution of the quadrupole force. To achieve for the 
quadrupole induced optical lattice the same depth as being produced for the dipole induced 
optical lattice potential in the NIST experiment one must increase the field intensity by 
4 orders of magnitude. This is not a practical intensity in nowadays technology. However, 
by increasing the intensity by less than 2 orders of magnitude, the depth of the quadrupole 
induced optical lattice potential is smaller just by about 2 orders of magnitude than in the 
NIST experiments. Such a potential can trap Na atoms provided that they are cooled down 
to less than 0.1/iK. The required temperature is sufficiently low to keep the atoms trapped 
in the quadrupole induced optical lattice potential well since /c#T < V max . 

Let us now turn to the case of Ca atoms where the first excited state is a D state. The 
results of our calculations which are presented in Fig. 2 show that when the laser light is 
detuned slightly from the S — > D transition (A = 457 nm, A^ = 70 kHz) it can trap the 
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Ca atoms solely by quadrupole interactions. The Ca atoms must, however, be pre-cooled. 
It is known that already today ultracold calcium can be 

this temperature the Ca atoms can be trapped by the quadrupole lattice potential using 
a laser intensity of 10 W/cm 2 . Other values of interest: In order to form a quadrupole 
induced optical lattice potential for field intensities which are smaller than 100 W/cm 2 , the 
Ca atoms should be cooled to a temperature which is smaller than about 30fiK. For the 
intensity of / = 1 W/cm 2 the Ca atoms will be trapped by the quadrupole force when 
their temperature is smaller than 3fiK. In passing we note that the used value of detuning 
Al = 70 kHz is relatively small compared to that of Fig. 1, but still remains more than 
200 times larger than the natural linewidth T « 300 Hz of the upper D state of Ca Q|. 
Therefore, the spontaneous emission effects can be safely neglected here. This statement 
has been confirmed also by an additional numerical calculation Q]. 

Above illustrative examples of course do not restrict the range of possible applications just 
to alkali metal or alkali earth atoms. Similar kind of quadrupole trapping can be exploited 
e.g. in the case of rare gases, which have been laser cooled and Bose condensed as well (see 

n 

for example Ref. |lfi|). 

In conclusion, we would like to re-stress that our formulation based upon the Floquet 
eigenstates is exact for all field intensities. This method provides a playground for a rational 
design of optical lattice potentials induced by the dipole or the quadrupole forces or by any 
combination of them. Design of quadrupole optical lattices might prove to be important in 
every practical situation where a specific lattice constant is required and where no dipole 
transition possesses a suitable wavelength. Such applications may include e.g. interference 
experiments with wavepackets of cold atoms, or an enhancement of the high harmonic gen- 
eration efficiency via constructive interference/phase matching. Besides trapping of cold 
atoms in optical lattices, the theoretical framework described in the present paper can be 
used also for other purposes, e.g. for studying an effect of quadrupole interaction on the 
performance of atomic clocks (cf. Ref. \^\)- 
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FIG. 1: The depth V mauX of the optical lattice potential well for Na atoms plotted as a function 
of the laser intensity /. The laser parameters, A = 343 nm and = 1 GHz, correspond to the 
transition from the S ground state to the lowest D state of Na. The solid line represents the 
numerically exact solutions of Eq. (J5J) which coalesces to all significant digits with the second order 
perturbation theory contribution of the quadrupole forces, see Eq. (J7J). The dashed line shows the 
contribution of dipole forces calculated from Eq. ©. 
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FIG. 2: The depth Vn iax of the optical lattice potential well for Ca atoms plotted as a function 
of the laser intensity /. The laser parameters, A = 457 nm and = 70 kHz, correspond to the 
transition from the S ground state to the first excited D state of Ca. The solid line stands here for 
the numerically exact optical lattice potential as obtained from the numerical solution of Eq. (J5J). 
The dashed line shows the contribution of quadrupole forces accounted within the second order 
perturbation theory, see Eq. (0). 
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